Recombinant SARS-CoV-2 lacking initiating and internal methionine codons within ORF10 is attenuated in vivo (preprint)

SARS-CoV-2 has been proposed to encode ORF10 as the 3' terminal gene in the viral genome. However, the potential role and even existence of a functional ORF10 product has been the subject of debate. There are significant structural features in the viral genomic RNA that could, by themselves, explain the retention of the ORF10 nucleotide sequences without the need for a functional protein product. To explore this question further we made two recombinant viruses, firstly a control virus (WT) based on the genome sequence of the original Wuhan isolate and with the inclusion of the early D614G mutation in the Spike protein. We also made a second virus, identical to WT except for two additional changes that replaced the initiating ORF10 start codon and an internal methionine codon for stop codons (ORF10KO). Here we show that the two viruses have apparently identical growth kinetics in a VeroE6 cell line that over expresses TMPRSS2 (VTN cells). However, in A549 cells over expressing ACE2 and TMPRSS2 (A549-AT cells) the ORF10KO virus appears to have a small growth rate advantage. Growth competition experiments were used whereby the two viruses were mixed, passaged in either VTN or A549-AT cells and the resulting output virus was sequenced. We found that in VTN cells the WT virus quickly dominated whereas in the A549-AT cells the ORF10KO virus dominated. We then used a hamster model of SARS-CoV-2 infection and determined that the ORF10KO virus has attenuated pathogenicity (as measured by weight loss). We found an almost 10-fold reduction in viral titre in the lower respiratory tract for ORF10KO vs WT. In contrast, the WT and ORF10KO viruses had similar titres in the upper respiratory tract. Sequencing of viral RNA in the lungs of hamsters infected with ORF10KO virus revealed that this virus frequently reverts to WT. Our data suggests that the retention of a functional ORF10 sequence is highly desirable for SARS-CoV-2 infection of hamsters and affects the virus's ability to propagate in the lower respiratory tract.

SARS-CoV-2 NSP12 associates with the TRiC complex and the P323L substitution is a host adaption. (preprint)

SARS-CoV-2 emerged into the human population in late 2019 and human to human transmission has dominated the evolutionary landscape and driven the selection of different lineages. The first major change that resulted in increased transmission was the D614G substitution in the spike protein. This was accompanied by the P323L substitution in the viral RNA dependent RNA polymerase (RdRp) (NSP12). Together, with D614G these changes are the root of the predominant global SARS-CoV-2 landscape. Here, we found that NSP12 formed an interactome with cellular proteins. The functioning of NSP12 was dependent on the T-complex protein Ring Complex, a molecular chaperone. In contrast, there was differential association between NSP12 variants and components of a phosphatase complex (PP2/PP2A and STRN3). Virus expressing NSP12L323 was less sensitive to perturbations in PP2A and supports the paradigm that ongoing genotype to phenotype adaptation of SARS-CoV-2 in humans is not exclusively restricted to the spike protein.

Development of SARS-CoV-2 replicons for the ancestral virus and variant of concern Delta for antiviral screening. (preprint)

SARS-CoV-2 is the aetiologic agent of COVID-19 and the associated ongoing pandemic. As the pandemic has progressed, Variants of Concern (VOC) have emerged with lineage defining mutations. Using a SARS-CoV-2 reverse genetic system, based on transformation associated recombination in yeast, a series of replicons were produced for the ancestral Wuhan virus and the SARS-CoV-2 VOC Delta in which different combinations of the Spike, membrane, ORF6 and ORF7a coding sequences were replaced with sequences encoding the selectable marker puromycin N-acetyl transferase and reporter proteins (Renilla luciferase, mNeonGreen and mScarlet). Replicon RNAs were replication competent in African green monkey kidney (Vero E6) derived cells and a range of human cell lines, with a Vero E6 cell line expressing ACE2 and TMPRSS2 showing much higher transfection efficiency and overall levels of Renilla luciferase activity. The replicons could be used for transient gene expression studies, but cell populations that stably maintained the replicons could not be propagated. Replication of the transiently expressed replicon RNA genomes was sensitive to remedesivir, providing a system to dissect the mechanism of action of antiviral compounds.

Impact of voluntary risk-mitigation behaviour on transmission of the Omicron SARS-CoV-2 variant in England (preprint)

Abstract Background The Omicron variant of SARS-CoV-2 infection poses substantial challenges to public health. In England, "plan B" mitigation measures were introduced in December 2021 including increased home working and face coverings in shops, but stopped short of restrictions on social contacts. The impact of voluntary risk mitigation behaviours on future SARS-CoV-2 burden is unknown. Methods We developed a rapid online survey of risk mitigation behaviours during the winter 2021 festive period and deployed in two longitudinal cohort studies in the UK (Avon Longitudinal Study of Parents and Children (ALSPAC) and TwinsUK/Covid Symptom Study (CSS) Biobank) in December 2021. Using an individual-based, probabilistic model of COVID-19 transmission between social contacts with SARS-CoV-2 Omicron variant parameters and realistic vaccine coverage in England, we describe the potential impact of the SARS-CoV-2 Omicron wave in England in terms of the effective reproduction number and cumulative infections, hospital admissions and deaths. Using survey results, we estimated in real-time the impact of voluntary risk mitigation behaviours on the Omicron wave in England, if implemented for the entire epidemic wave. Results Over 95% of survey respondents (N_ALSPAC=2,686 and N_Twins=6,155) reported some risk mitigation behaviours, with being fully vaccinated and using home testing kits the most frequently reported behaviours. Less than half of those respondents reported that their behaviour was due to "plan B". We estimate that without risk mitigation behaviours, the Omicron variant is consistent with an effective reproduction number between 2.5 and 3.5. Due to the reduced vaccine effectiveness against infection with the Omicron variant, our modelled estimates suggest that between 55% and 60% of the English population could be infected during the current wave, translating into between 15,000 and 46,000 cumulative deaths, depending on assumptions about vaccine effectiveness. We estimate that voluntary risk reduction measures could reduce the effective reproduction number to between 1.8 and 2.2 and reduce the cumulative number of deaths by up to 24%. Conclusions We conclude that voluntary measures substantially reduce the projected impact of the SARS-CoV-2 Omicron variant, but that voluntary measures alone would be unlikely to completely control transmission.

Rapid selection of P323L in the SARS-CoV-2 polymerase (NSP12) in humans and non-human primate models and confers a large plaque phenotype (preprint)

The mutational landscape of SARS-CoV-2 varies at both the dominant viral genome sequence and minor genomic variant population. An early change associated with transmissibility was the D614G substitution in the spike protein. This appeared to be accompanied by a P323L substitution in the viral polymerase (NSP12), but this latter change was not under strong selective pressure. Investigation of P323L/D614G changes in the human population showed rapid emergence during the containment phase and early surge phase of wave 1 in the UK. This rapid substitution was from minor genomic variants to become part of the dominant viral genome sequence. A rapid emergence of 323L but not 614G was observed in a non-human primate model of COVID-19 using a starting virus with P323 and D614 in the dominant genome sequence and 323L and 614G in the minor variant population. In cell culture, a recombinant virus with 323L in NSP12 had a larger plaque size than the same recombinant virus with P323. These data suggest that it may be possible to predict the emergence of a new variant based on tracking the distribution and frequency of minor variant genomes at a population level, rather than just focusing on providing information on the dominant viral genome sequence e.g., consensus level reporting. The ability to predict an emerging variant of SARS-CoV-2 in the global landscape may aid in the evaluation of medical countermeasures and non-pharmaceutical interventions.

Nanopore ReCappable Sequencing maps SARS-CoV-2 5' capping sites and provides new insights into the structure of sgRNAs (preprint)

The SARS-CoV-2 virus has a complex transcriptome characterised by multiple, nested sub genomic RNAs used to express structural and accessory proteins. Long-read sequencing technologies such as nanopore direct RNA sequencing can recover full-length transcripts, greatly simplifying the assembly of structurally complex RNAs. However, these techniques do not detect the 5' cap, thus preventing reliable identification and quantification of full-length, coding transcript models. Here we used Nanopore ReCappable Sequencing (NRCeq), a new technique that can identify capped full-length RNAs, to assemble a complete annotation of SARS-CoV-2 sgRNAs and annotate the location of capping sites across the viral genome. We obtained robust estimates of sgRNA expression across cell lines and viral isolates and identified novel canonical and non-canonical sgRNAs, including one that uses a previously un-annotated leader-to-body junction site. The data generated in this work constitute a useful resource for the scientific community and provide important insights into the mechanisms that regulate the transcription of SARS-CoV-2 sgRNAs.

Chronic SARS-CoV-2 infection and viral evolution in a hypogammaglobulinaemic individual. (preprint)

There is widespread interest in the capacity for SARS-CoV-2 evolution in the face of selective pressures from host immunity, either naturally acquired post-exposure or from vaccine acquired immunity. Allied to this is the potential for long perm persistent infections within immune compromised individuals to allow a broader range of viral evolution in the face of sub-optimal immune driven selective pressure. Here we report on an immunocompromised individual who is hypogammaglobulinaemic and was persistently infected with SARS-CoV-2 for over 290 days, the longest persistent infection recorded in the literature to date. During this time, nine samples of viral nucleic acid were obtained and analysed by next-generation sequencing. Initially only a single mutation (L179I) was detected in the spike protein relative to the prototypic SARS-CoV-2 Wuhan-Hu-1 isolate, with no further changes identified at day 58. However, by day 155 the spike protein had acquired a further four amino acid changes, namely S255F, S477N, H655Y and D1620A and a two amino acid deletion ({Delta}H69/{Delta}V70). Infectious virus was cultured from a nasopharyngeal sample taken on day 155 and next-generation sequencing confirmed that the mutations in the virus mirrored those identified by sequencing of the corresponding swab sample. The isolated virus was susceptible to remdesivir in vitro, however a 17-day course of remdesivir started on day 213 had no effect on the viral RT-PCR cycle threshold (Ct) value. On day 265 the patient was treated with the combination of casirivimab and imdevimab. The patient experienced progressive resolution of all symptoms over the next 8 weeks and by day 311 the virus was no longer detectable by RT-PCR. The {Delta}H69/{Delta}V70 deletion in the N-terminus of the spike protein which arose in our patient is also present in the B.1.1.7 variant of concern and has been associated with viral escape mutagenesis after treatment of another immunocompromised patient with convalescent plasma. Our data confirms the significance of this deletion in immunocompromised patients but illustrates it can arise independently of passive antibody transfer, suggesting the deletion may be an enabling mutation that compensates for distant changes in the spike protein that arise under selective pressure.

Structure and mechanism of SARS-CoV-2 Spike N679-V687 deletion variant elucidate cell-type specific evolution of viral fitness (preprint)

As the global burden of SARS-CoV-2 infections escalates, so does the evolution of viral variants which is of particular concern due to their potential for increased transmissibility and pathology. In addition to this entrenched variant diversity in circulation, RNA viruses can also display genetic diversity within single infected hosts with co-existing viral variants evolving differently in distinct cell types. The BriS{Delta} variant, originally identified as a viral subpopulation by passaging SARS-CoV-2 isolate hCoV-19/England/02/2020, comprises in the spike glycoprotein an eight amino-acid deletion encompassing the furin recognition motif and S1/S2 cleavage site. Here, we analyzed the structure, function and molecular dynamics of this variant spike, providing mechanistic insight into how the deletion correlates to viral cell tropism, ACE2 receptor binding and infectivity, allowing the virus to probe diverse trajectories in distinct cell types to evolve viral fitness. TeaserSARS-CoV-2 can exploit different cell types to diversify and evolve virus variants distinct in infectivity and structure.

Identification and quantification of SARS-CoV-2 leader subgenomic mRNA gene junctions in nasopharyngeal samples shows phasic transcription in animal models of COVID-19 and aberrant pattens in humans (preprint)

Introduction: SARS-CoV-2 has a complex strategy for the transcription of viral subgenomic mRNAs (sgmRNAs), which are targets for nucleic acid diagnostics. Each of these sgRNAs has a unique 5 sequence, the leader-transcriptional regulatory sequence gene junction (leader-TRS-junction), that can be identified using sequencing. Results: High resolution sequencing has been used to investigate the biology of SARS-CoV-2 and the host response in cell culture models and from clinical samples. LeTRS, a bioinformatics tool, was developed to identify leader-TRS-junctions and be used as a proxy to quantify sgmRNAs for understanding virus biology. This was tested on published datasets and clinical samples from patients and longitudinal samples from animal models with COVID-19. Discussion: LeTRS identified known leader-TRS-junctions and identified novel species that were common across different species. The data indicated multi-phasic abundance of sgmRNAs in two different animal models, with spikes in sgmRNA abundance reflected in human samples, and therefore has implications for transmission models and nucleic acid-based diagnostics.

SARS-CoV-2 candidate vaccine ChAdOx1 nCoV-19 infection of human cell lines reveals a normal low range of viral backbone gene expression alongside very high levels of SARS-CoV-2 S glycoprotein expression. (preprint)

Background: ChAdOx1 nCoV-19 is a recombinant adenovirus vaccine candidate against SARS-CoV-2. Although replication defective in normal cells, 28kbp of adenovirus genes are delivered to the cell nucleus alongside the SARS-CoV-2 S glycoprotein gene.Methods: We used direct RNA sequencing to analyse transcript expression from the ChAdOx1 nCoV-19 genome in human MRC-5 and A549 cell lines that are non-permissive for vector replication alongside the replication permissive cell line, HEK293. In addition, we used quantitative proteomics to study over time the proteome and phosphoproteome of A549 and MRC5 cells infected with the ChAdOx1 nCoV-19 vaccine candidate.Results: The expected SARS-CoV-2 S coding transcript dominated in all cell lines. We also detected rare S transcripts with aberrant splice patterns or polyadenylation site usage. Adenovirus vector transcripts were almost absent in MRC-5 cells but in A549 cells there was a broader repertoire of adenoviral gene expression at very low levels. Proteomically, in addition to S glycoprotein, we detected multiple adenovirus proteins in A549 cells compared to just one in MRC5 cells. Conclusions: Overall the ChAdOx1 nCoV-19 vaccine’s transcriptomic and proteomic repertoire is as expected. The combined transcriptomic and proteomics approaches provide an unparalleled insight into the behaviour of this important class of vaccine candidate and illustrate the potential of this technique to inform future viral vaccine vector design.

Sequential infection with influenza A virus followed by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) leads to more severe disease and encephalitis in a mouse model of COVID-19. (preprint)

COVID-19 is a spectrum of clinical symptoms in humans caused by infection with SARS-CoV-2, a recently emerged coronavirus that has rapidly caused a pandemic. Coalescence of a second wave of this virus with seasonal respiratory viruses, particularly influenza virus is a possible global health concern. To investigate this, transgenic mice expressing the human ACE2 receptor driven by the epithelial cell cytokeratin-18 gene promoter (K18-hACE2) were first infected with IAV followed by SARS-CoV-2. The host response and effect on virus biology was compared to K18-hACE2 mice infected with IAV or SARS-CoV-2 only. Infection of mice with each individual virus resulted in a disease phenotype compared to control mice. Although, SARS-CoV-2 RNA synthesis appeared significantly reduced in the sequentially infected mice, these mice had a more rapid weight loss, more severe lung damage and a prolongation of the innate response compared to singly infected or control mice. The sequential infection also exacerbated the extrapulmonary manifestations associated with SARS-CoV-2. This included a more severe encephalitis. Taken together, the data suggest that the concept of "twinfection" is deleterious and mitigation steps should be instituted as part of a comprehensive public health response to the COVID-19 pandemic.

Neuropilin-1 is a host factor for SARS-CoV-2 infection (preprint)

SARS-CoV-2 is the causative agent of COVID-19, a coronavirus disease that has infected more than 6.6 million people and caused over 390,000 deaths worldwide1,2. The Spike (S) protein of the virus forms projections on the virion surface responsible for host cell attachment and penetration. This viral glycoprotein is synthesized as a precursor in infected cells and, to be active, must be cleaved to two associated polypeptides: S1 and S2(3,4). For SARS-CoV-2 the cleavage is catalysed by furin, a host cell protease, which cleaves the S protein precursor at a specific sequence motif that generates a polybasic Arg-Arg-Ala-Arg (RRAR) C-terminal sequence on S1. This sequence motif conforms to the C-end rule (CendR), which means that the C-terminal sequence may allow the protein to associate with cell surface neuropilin-1 (NRP1) and neuropilin-2 (NRP2) receptors5. Here we demonstrate using immunoprecipitation, site-specific mutagenesis, structural modelling, and antibody blockade that, in addition to engaging the known receptor ACE2, S1 can bind to NRP1 through the canonical CendR mechanism. This interaction enhances infection by SARS-CoV-2 in cell culture. NRP1 thus serves as a host factor for SARS-CoV-2 infection, and provides a therapeutic target for COVID-19.

Characterisation of the transcriptome and proteome of SARS-CoV-2 using direct RNA sequencing and tandem mass spectrometry reveals evidence for a cell passage induced in-frame deletion in the spike glycoprotein that removes the furin-like cleavage site. (preprint)

Direct RNA sequencing using an Oxford Nanopore MinION characterised the transcriptome of SARS-CoV-2 grown in Vero E6 cells. This cell line is being widely used to propagate the novel coronavirus. The viral transcriptome was analysed using a recently developed ORF-centric pipeline. This revealed the pattern of viral transcripts, (i.e. subgenomic mRNAs), generally fitted the predicted replication and transcription model for coronaviruses. A 24 nt in-frame deletion was detected in subgenomic mRNAs encoding the spike (S) glycoprotein. This feature was identified in over half of the mapped transcripts and was predicted to remove a proposed furin cleavage site from the S glycoprotein. This motif directs cleavage of the S glycoprotein into functional subunits during virus entry or exit. Cleavage of the S glycoprotein can be a barrier to zoonotic coronavirus transmission and affect viral pathogenicity. Allied to this transcriptome analysis, tandem mass spectrometry was used to identify over 500 viral peptides and 44 phosphopeptides, covering almost all of the proteins predicted to be encoded by the SARS-CoV-2 genome, including peptides unique to the deleted variant of the S glycoprotein. Detection of an apparently viable deletion in the furin cleavage site of the S glycoprotein reinforces the point that this and other regions of SARS-CoV-2 proteins may readily mutate. This is of clear significance given the interest in the S glycoprotein as a potential vaccine target and the observation that the furin cleavage site likely contributes strongly to the pathogenesis and zoonosis of this virus. The viral genome sequence should be carefully monitored during the growth of viral stocks for research, animal challenge models and, potentially, in clinical samples. Such variations may result in different levels of virulence, morbidity and mortality.

Amplicon based MinION sequencing of SARS-CoV-2 and metagenomic characterisation of nasopharyngeal swabs from patients with COVID-19 (preprint)

COVID-19 is a complex disease phenotype where the underlying microbiome could influence morbidity and mortality. Amplicon and metagenomic MinION based sequencing was used to rapidly (within 8 hours) identify SARS-CoV-2 and assess the microbiome in nasopharyngeal swabs obtained from patients with COVID-19 by the ISARIC 4C consortium.